Simultaneous profiling of chromatin accessibility and DNA methylation in complete plant genomes using long-read sequencing.

Epigenetic regulations, including chromatin accessibility, nucleosome positioning and DNA methylation intricately shape genome function. However, current chromatin profiling techniques relying on short-read sequencing technologies fail to characterise highly repetitive genomic regions and cannot detect multiple chromatin features simultaneously. Here, we performed Simultaneous Accessibility and DNA Methylation Sequencing (SAM-seq) of purified plant nuclei. Thanks to the use of long-read nanopore sequencing, SAM-seq enables high-resolution profiling of m6A-tagged chromatin accessibility together with endogenous cytosine methylation in plants. Analysis of naked genomic DNA revealed significant sequence preference biases of m6A-MTases, controllable through a normalisation step. By applying SAM-seq to Arabidopsis and maize nuclei we obtained fine-grained accessibility and DNA methylation landscapes genome-wide. We uncovered crosstalk between chromatin accessibility and DNA methylation within nucleosomes of genes, TEs, and centromeric repeats. SAM-seq also detects DNA footprints over cis-regulatory regions. Furthermore, using the single-molecule information provided by SAM-seq we identified extensive cellular heterogeneity at chromatin domains with antagonistic chromatin marks, suggesting that bivalency reflects cell-specific regulations. SAM-seq is a powerful approach to simultaneously study multiple epigenetic features over unique and repetitive sequences, opening new opportunities for the investigation of epigenetic mechanisms.

Pervasive G × E interactions shape adaptive trajectories and the exploration of the phenotypic space in artificial selection experiments.

Quantitative genetics models have shown that long-term selection responses depend on initial variance and mutational influx. Understanding limits of selection requires quantifying the role of mutational variance. However, correlative responses to selection on nonfocal traits can perturb the selection response on the focal trait; and generations are often confounded with selection environments so that genotype by environment (G×E) interactions are ignored. The Saclay divergent selection experiments (DSEs) on maize flowering time were used to track the fate of individual mutations combining genotyping data and phenotyping data from yearly measurements (DSEYM) and common garden experiments (DSECG) with four objectives: (1) to quantify the relative contribution of standing and mutational variance to the selection response, (2) to estimate genotypic mutation effects, (3) to study the impact of G×E interactions in the selection response, and (4) to analyze how trait correlations modulate the exploration of the phenotypic space. We validated experimentally the expected enrichment of fixed beneficial mutations with an average effect of +0.278 and +0.299 days to flowering, depending on the genetic background. Fixation of unfavorable mutations reached up to 25% of incoming mutations, a genetic load possibly due to antagonistic pleiotropy, whereby mutations fixed in the selection environment (DSEYM) turned to be unfavorable in the evaluation environment (DSECG). Global patterns of trait correlations were conserved across genetic backgrounds but exhibited temporal patterns. Traits weakly or uncorrelated with flowering time triggered stochastic exploration of the phenotypic space, owing to microenvironment-specific fixation of standing variants and pleiotropic mutational input.

Identification of Key Tissue-Specific, Biological Processes by Integrating Enhancer Information in Maize Gene Regulatory Networks.

Enhancers are key players in the spatio-temporal coordination of gene expression during numerous crucial processes, including tissue differentiation across development. Characterizing the transcription factors (TFs) and genes they connect, and the molecular functions underpinned is important to better characterize developmental processes. In plants, the recent molecular characterization of enhancers revealed their capacity to activate the expression of several target genes. Nevertheless, identifying these target genes at a genome-wide level is challenging, particularly for large-genome species, where enhancers and target genes can be hundreds of kilobases away. Therefore, the contribution of enhancers to plant regulatory networks remains poorly understood. Here, we investigate the enhancer-driven regulatory network of two maize tissues at different stages: leaves at seedling stage (V2-IST) and husks (bracts) at flowering. Using systems biology, we integrate genomic, epigenomic, and transcriptomic data to model the regulatory relationships between TFs and their potential target genes, and identify regulatory modules specific to husk and V2-IST. We show that leaves at the V2-IST stage are characterized by the response to hormones and macromolecules biogenesis and assembly, which are regulated by the BBR/BPC and AP2/ERF TF families, respectively. In contrast, husks are characterized by cell wall modification and response to abiotic stresses, which are, respectively, orchestrated by the C2C2/DOF and AP2/EREB families. Analysis of the corresponding enhancer sequences reveals that two different transposable element families (TIR transposon Mutator and MITE Pif/Harbinger) have shaped part of the regulatory network in each tissue, and that MITEs have provided potential new TF binding sites involved in husk tissue-specificity.

High throughput genotyping of structural variations in a complex plant genome using an original Affymetrix® axiom® array.

BACKGROUND: Insertions/deletions (InDels) and more specifically presence/absence variations (PAVs) are pervasive in several species and have strong functional and phenotypic effect by removing or drastically modifying genes. Genotyping of such variants on large panels remains poorly addressed, while necessary for approaches such as association mapping or genomic selection. RESULTS: We have developed, as a proof of concept, a new high-throughput and affordable approach to genotype InDels. We first identified 141,000 InDels by aligning reads from the B73 line against the genome of three temperate maize inbred lines (F2, PH207, and C103) and reciprocally. Next, we designed an Affymetrix® Axiom® array to target these InDels, with a combination of probes selected at breakpoint sites (13%) or within the InDel sequence, either at polymorphic (25%) or non-polymorphic sites (63%) sites. The final array design is composed of 662,772 probes and targets 105,927 InDels, including PAVs ranging from 35 bp to 129kbp. After Affymetrix® quality control, we successfully genotyped 86,648 polymorphic InDels (82% of all InDels interrogated by the array) on 445 maize DNA samples with 422,369 probes. Genotyping InDels using this approach produced a highly reliable dataset, with low genotyping error (~ 3%), high call rate (~ 98%), and high reproducibility (> 95%). This reliability can be further increased by combining genotyping of several probes calling the same InDels (< 0.1% error rate and > 99.9% of call rate for 5 probes). This "proof of concept" tool was used to estimate the kinship matrix between 362 maize lines with 57,824 polymorphic InDels. This InDels kinship matrix was highly correlated with kinship estimated using SNPs from Illumina 50 K SNP arrays. CONCLUSIONS: We efficiently genotyped thousands of small to large InDels on a sizeable number of individuals using a new Affymetrix® Axiom® array. This powerful approach opens the way to studying the contribution of InDels to trait variation and heterosis in maize. The approach is easily extendable to other species and should contribute to decipher the biological impact of InDels at a larger scale.

Sequence analysis of European maize inbred line F2 provides new insights into molecular and chromosomal characteristics of presence/absence variants.

BACKGROUND: Maize is well known for its exceptional structural diversity, including copy number variants (CNVs) and presence/absence variants (PAVs), and there is growing evidence for the role of structural variation in maize adaptation. While PAVs have been described in this important crop species, they have been only scarcely characterized at the sequence level and the extent of presence/absence variation and relative chromosomal landscape of inbred-specific regions remain to be elucidated. RESULTS: De novo genome sequencing of the French F2 maize inbred line revealed 10,044 novel genomic regions larger than 1 kb, making up 88 Mb of DNA, that are present in F2 but not in B73 (PAV). This set of maize PAV sequences allowed us to annotate PAV content and to analyze sequence breakpoints. Using PAV genotyping on a collection of 25 temperate lines, we also analyzed Linkage Disequilibrium in PAVs and flanking regions, and PAV frequencies within maize genetic groups. CONCLUSIONS: We highlight the possible role of MMEJ-type double strand break repair in maize PAV formation and discover 395 new genes with transcriptional support. Pattern of linkage disequilibrium within PAVs strikingly differs from this of flanking regions and is in accordance with the intuition that PAVs may recombine less than other genomic regions. We show that most PAVs are ancient, while some are found only in European Flint material, thus pinpointing structural features that may be at the origin of adaptive traits involved in the success of this material. Characterization of such PAVs will provide useful material for further association genetic studies in European and temperate maize.

Independent introductions and admixtures have contributed to adaptation of European maize and its American counterparts.

Through the local selection of landraces, humans have guided the adaptation of crops to a vast range of climatic and ecological conditions. This is particularly true of maize, which was domesticated in a restricted area of Mexico but now displays one of the broadest cultivated ranges worldwide. Here, we sequenced 67 genomes with an average sequencing depth of 18x to document routes of introduction, admixture and selective history of European maize and its American counterparts. To avoid the confounding effects of recent breeding, we targeted germplasm (lines) directly derived from landraces. Among our lines, we discovered 22,294,769 SNPs and between 0.9% to 4.1% residual heterozygosity. Using a segmentation method, we identified 6,978 segments of unexpectedly high rate of heterozygosity. These segments point to genes potentially involved in inbreeding depression, and to a lesser extent to the presence of structural variants. Genetic structuring and inferences of historical splits revealed 5 genetic groups and two independent European introductions, with modest bottleneck signatures. Our results further revealed admixtures between distinct sources that have contributed to the establishment of 3 groups at intermediate latitudes in North America and Europe. We combined differentiation- and diversity-based statistics to identify both genes and gene networks displaying strong signals of selection. These include genes/gene networks involved in flowering time, drought and cold tolerance, plant defense and starch properties. Overall, our results provide novel insights into the evolutionary history of European maize and highlight a major role of admixture in environmental adaptation, paralleling recent findings in humans.

The bright side of transposons in crop evolution.

The past decades have revealed an unexpected yet prominent role of so-called 'junk DNA' in the regulation of gene expression, thereby challenging our view of the mechanisms underlying phenotypic evolution. In particular, several mechanisms through which transposable elements (TEs) participate in functional genome diversity have been depicted, bringing to light the 'TEs bright side'. However, the relative contribution of those mechanisms and, more generally, the importance of TE-based polymorphisms on past and present phenotypic variation in crops species remain poorly understood. Here, we review current knowledge on both issues, and discuss how analyses of massively parallel sequencing data combined with statistical methodologies and functional validations will help unravelling the impact of TEs on crop evolution in a near future.

Transposable elements, a treasure trove to decipher epigenetic variation: insights from Arabidopsis and crop epigenomes.

In the past decade, plant biologists and breeders have developed a growing interest in the field of epigenetics, which is defined as the study of heritable changes in gene expression that cannot be explained by changes in the DNA sequence. Epigenetic marks can be responsive to the environment, and evolve faster than genetic changes. Therefore, epigenetic diversity may represent an unexplored resource of natural variation that could be used in plant breeding programmes. On the other hand, crop genomes are largely populated with transposable elements (TEs) that are efficiently targeted by epigenetic marks, and part of the epigenetic diversity observed might be explained by TE polymorphisms. Characterizing the degree to which TEs influence epigenetic variation in crops is therefore a major goal to better use epigenetic variation. To date, epigenetic analyses have been mainly focused on the model plant Arabidopsis thaliana, and have provided clues on epigenome features, components that silence pathways, and effects of silencing impairment. But to what extent can Arabidopsis be used as a model for the epigenomics of crops? In this review, we discuss the similarities and differences between the epigenomes of Arabidopsis and crops. We explore the relationship between TEs and epigenomes, focusing on TE silencing control and escape, and the impact of TE mobility on epigenomic variation. Finally, we provide insights into challenges to tackle, and future directions to take in the route towards using epigenetic diversity in plant breeding programmes.

Vernalization treatment induces site-specific DNA hypermethylation at the VERNALIZATION-A1 (VRN-A1) locus in hexaploid winter wheat.

BACKGROUND: Certain temperate species require prolonged exposure to low temperature to initiate transition from vegetative growth to flowering, a process known as vernalization. In wheat, winter cultivars require vernalization to initiate flowering, making vernalization requirement a trait of key importance in wheat agronomy. The genetic bases of vernalization response have been largely studied in wheat, leading to the characterization of a regulation pathway that involves the key gene VERNALIZATION1 (VRN1). While previous studies in wheat and barley have revealed the functional role of histone modification in setting VRN1 expression, other mechanisms might also be involved. Here, we were interested in determining whether the cold-induced expression of the wheat VRN-A1 gene is associated with a change in DNA methylation. RESULTS: We provide the first DNA methylation analysis of the VRN-A1 gene, and describe the existence of methylation at CG but also at non CG sites. While CG sites show a bell-shape profile typical of gene-body methylation, non CG methylation is restricted to the large (8.5 kb) intron 1, in a region harboring fragments of transposable elements (TEs). Interestingly, cold induces a site-specific hypermethylation at these non CG sites. This increase in DNA methylation is transmitted through mitosis, and is reset to its original level after sexual reproduction. CONCLUSIONS: These results demonstrate that VRN-A1 has a particular DNA methylation pattern, exhibiting rapid shift within the life cycle of a winter wheat plant following exposure to particular environmental conditions. The finding that this shift occurs at non CG sites in a TE-rich region opens interesting questions onto the possible consequences of this type of methylation in gene expression.

International Congress on Transposable Elements (ICTE) 2012 in Saint Malo and the sea of TE stories.

An international conference on Transposable Elements (TEs) was held 21-24 April 2012 in Saint Malo, France. Organized by the French Transposition Community (GDR Elements Génétiques Mobiles et Génomes, CNRS) and the French Society of Genetics (SFG), the conference's goal was to bring together researchers from around the world who study transposition in diverse organisms using multiple experimental approaches. The meeting drew more than 217 attendees and most contributed through poster presentations (117), invited talks and short talks selected from poster abstracts (48 in total). The talks were organized into four scientific sessions, focused on: impact of TEs on genomes, control of transposition, evolution of TEs and mechanisms of transposition. Here, we present highlights from the talks given during the platform sessions. The conference was sponsored by Alliance pour les sciences de la vie et de la santé (Aviesan), Centre national de la recherche scientifique (CNRS), Institut national de la santé et de la recherche médicale (INSERM), Institut de recherche pour le développement (IRD), Institut national de la recherche agronomique (INRA), Université de Perpignan, Université de Rennes 1, Région Bretagne and Mobile DNA. CHAIR OF THE ORGANIZATION COMMITTEE: Jean-Marc Deragon ORGANIZERS: Abdelkader Ainouche, Mireille Bétermier, Mick Chandler, Richard Cordaux, Gaël Cristofari, Jean-Marc Deragon, Pascale Lesage, Didier Mazel, Olivier Panaud, Hadi Quesneville, Chantal Vaury, Cristina Vieira and Clémentine Vitte.

A GeneTrek analysis of the maize genome.

Analysis of the sequences of 74 randomly selected BACs demonstrated that the maize nuclear genome contains approximately 37,000 candidate genes with homologues in other plant species. An additional approximately 5,500 predicted genes are severely truncated and probably pseudogenes. The distribution of genes is uneven, with approximately 30% of BACs containing no genes. BAC gene density varies from 0 to 7.9 per 100 kb, whereas most gene islands contain only one gene. The average number of genes per gene island is 1.7. Only 72% of these genes show collinearity with the rice genome. Particular LTR retrotransposon families (e.g., Gyma) are enriched on gene-free BACs, most of which do not come from pericentromeres or other large heterochromatic regions. Gene-containing BACs are relatively enriched in different families of LTR retrotransposons (e.g., Ji). Two major bursts of LTR retrotransposon activity in the last 2 million years are responsible for the large size of the maize genome, but only the more recent of these is well represented in gene-containing BACs, suggesting that LTR retrotransposons are more efficiently removed in these domains. The results demonstrate that sample sequencing and careful annotation of a few randomly selected BACs can provide a robust description of a complex plant genome.

LTR retrotransposons in rice (Oryza sativa, L.): recent burst amplifications followed by rapid DNA loss.

BACKGROUND: LTR retrotransposons are one of the main causes for plant genome size and structure evolution, along with polyploidy. The characterization of their amplification and subsequent elimination of the genomes is therefore a major goal in plant evolutionary genomics. To address the extent and timing of these forces, we performed a detailed analysis of 41 LTR retrotransposon families in rice. RESULTS: Using a new method to estimate the insertion date of both truncated and complete copies, we estimated these two forces more accurately than previous studies based on other methods. We show that LTR retrotransposons have undergone bursts of amplification within the past 5 My. These bursts vary both in date and copy number among families, revealing that each family has a particular amplification history. The number of solo LTR varies among families and seems to correlate with LTR size, suggesting that solo LTR formation is a family-dependent process. The deletion rate estimate leads to the prediction that the half-life of LTR retrotransposon sequences evolving neutrally is about 19 My in rice, suggesting that other processes than the formation of small deletions are prevalent in rice DNA removal. CONCLUSION: Our work provides insights into the dynamics of LTR retrotransposons in the rice genome. We show that transposable element families have distinct amplification patterns, and that the turn-over of LTR retrotransposons sequences is rapid in the rice genome.

Analysis of retrotransposon structural diversity uncovers properties and propensities in angiosperm genome evolution.

Analysis of LTR retrotransposon structures in five diploid angiosperm genomes uncovered very different relative levels of different types of genomic diversity. All species exhibited recent LTR retrotransposon mobility and also high rates of DNA removal by unequal homologous recombination and illegitimate recombination. The larger plant genomes contained many LTR retrotransposon families with >10,000 copies per haploid genome, whereas the smaller genomes contained few or no LTR retrotransposon families with >1,000 copies, suggesting that this differential potential for retroelement amplification is a primary factor in angiosperm genome size variation. The average ratios of transition to transversion mutations (Ts/Tv) in diverging LTRs were >1.5 for each species studied, suggesting that these elements are mostly 5-methylated at cytosines in an epigenetically silenced state. However, the diploid wheat Triticum monococcum and barley have unusually low Ts/Tv values (respectively, 1.9 and 1.6) compared with maize (3.9), medicago (3.6), and lotus (2.5), suggesting that this silencing is less complete in the two Triticeae. Such characteristics as the ratios of point mutations to indels (insertions and deletions) and the relative efficiencies of DNA removal by unequal homologous recombination compared with illegitimate recombination were highly variable between species. These latter variations did not correlate with genome size or phylogenetic relatedness, indicating that they frequently change during the evolutionary descent of plant lineages. In sum, the results indicate that the different sizes, contents, and structures of angiosperm genomes are outcomes of the same suite of mechanistic processes, but acting with different relative efficiencies in different plant lineages.